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Abstract 


This work proposes the use of a rigorous approach to the analysis of the fuel-cell diffusive resistances not only at the commercial scale, but also 
at the laboratory one. The attention used experimentally for checking diffusion resistances in electrodes, cells and stacks should imply the same 


attention in the corresponding data analysis techniques. 


For this reason, some corrections to the most common interpretation and correlation procedures have been introduced on the basis of the study 


of complex electrodes and comparing these ones with simple ones. 


In Part 1 of this work a theoretical discussion of some simple reference cases will be reported. In Part 2 some examples of applications to the 


interpretation of the experimental data will be given. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


The experimental analysis of heterogeneous kinetics is based, 
on the one hand, on very sophisticated measurement methods 
while, on the other, it often uses simplified methods for the 
interpretation of the results. In particular, for electrochemical 
kinetics, this involves the analysis of the complex electrodes 
working on non-uniform concentration fields: the current- 
voltage experimental results are often just interpreted in the same 
manner as those of the simple electrodes, that is in terms of two 
parameters only, the exchange current and the limit current, and 
without any particular attention to the effects of the geometry, 
the flow pattern and the consequent distribution of concentration 
on the active surface. 

The complexity of cells and stacks at the commercial scale 
is well known: ever more detailed and analytical answers to 
the problem can be found in detailed simulation tools [1-7]. 
But the detailed models often use data from laboratory-scale 
measurements [8,9], whose interpretation can be affected by the 
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scarce consideration of the complex phenomenology that also 
exists at relatively small scales. For this reason, a brief review 
of the recent history of applied research on fuel cells could 
be useful: in research laboratories, and specifically industrial 
ones, attention is being shifted more and more frequently from 
centimetre cell sizes (“postage stamp” scale) to decimetre cell 
sizes (“postcard” scale). The traditional methods for analysing 
the electrochemical kinetics, which have been elaborated for 
the smaller scale and applied to it with considerable success, 
are often too simplified for the larger scale, where the unifor- 
mity hypothesis (simple electrodes) is very often inadequate. 
In brief, the transition from the postage stamp to the postcard 
scale requires some changes in the handling of the laboratory 
data. 

Similar situations, where the points of view can be articulated 
at a local level, or scale, an aggregate level and a macroscopic 
level, have been encountered in many different applicative fields, 
but perhaps the relationships between the levels have been better 
studied and understood for the heterogeneous chemical reactions 
of traditional chemical engineering (catalytic or non-catalytic 
fluid—solid reactions, fluid—fluid reactions, slurries, trickle beds, 
etc. [10]). From this wide scenario of cases, which are, for the 
major part, well established in the field of chemical engineering, 
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Nomenclature 


Noe eH SS CRNA RTA 


corrective coefficient, from the averaging opera- 
tion (—) 

see Eq. (B.5) (-) 

see Eq. (B.7) (-) 

volumetric concentration of the component i 
(kmol/m?) 

see Eq. (B.9) (-) 

diffusivity (m/s) 

see Eq. (B.4) (-) 

Faraday’s constant (A s/mol) 

current density (A/m?) 

mean electrode current density (A/m?) 

kinetic constants (depending on the reaction 
orders) 

apparent kinetic constant (m/s) 

mass transfer coefficient (m/s) 

flow dynamic coefficient (m/s) 

intrinsic kinetic constant (m/s) 

global kinetic constant, local value (m/s) 

length (m) 

stoichiometric coefficient for electrons, absolute 
value (—) 

volumetric flow rate (m3/s) 

reaction rate, per unit surface (kmol/m? s) 

limit reaction rate (kmol/m? s) 

mean reaction rate (kmol/m? s) 

gas constant (j/mol K) 

thickness (m) 

electrode surface (m?) 

absolute temperature (K) 

utilisation factor (—) 

velocity (m/s) 

(in Appendix A) adimensional voltage (—) 
voltage (W/A) 

reaction rate (kmol/m?) 

see Eq. (B.1) (-) 

spatial coordinate, orthogonal to the electrode (m) 
see Eq. (B.2) (-) 

spatial coordinate, parallel to the electrode (m) 


Greek symbols 


Qj 
Bi 
B 
n 
p 
(e 


Vi 


g 


reaction orders, direct reaction (—) 
reaction orders, inverse reaction (—) 
symmetry factor (—) 

a-dimensional voltage loss (—) 
effectiveness factor (—) 

Thiele modulus (—) 

stoichiometric coefficient (—) 
adimensional axial coordinate (—) 


Apexes and indexes 


A 


key reagent (hydrogen at the anode, oxygen at the 
cathode) 
concentration 


eff effective 

eq equilibrium 

ex external 

i generic component 
L limit 

m mean value 

o inlet 

r reference 

s electrode surface 

* exchange 

# refers to the shifting from the diagonal of Eq. (B.4) 
£ 


refers to the other electrode 


some useful approaches can be drawn for an analogous study of 
many cases of electrochemical interest. 

In analysing complex electrodes, overcoming any immedi- 
ate and substantially semi-empirical and superficial approach 
requires an advanced study of the transport phenomena of 
momentum, mass, heat and charge and the real electrode geom- 
etry and structure, so that non-trivial physical—-mathematical 
simulation approaches and their related numerical tools are nec- 
essary at the laboratory scale [11]. In other words, a wider point 
of view is needed, through which, on the one hand, the level 
of complexity necessary and sufficient for a good analysis of 
a set of experimental data can be determined and, on the other 
hand, the correct formal analytical tools, at the proper level of 
simplification, can be found. 

Of particular interest, and particular difficulty too, is an anal- 
ysis that highlights the behaviour of a single electrode of a couple 
of complex electrodes. In such cases it is necessary to maintain 
the behaviour of the second electrode constant and uniform, 
approximating the behaviour of a simple electrode, so that it can 
be considered non-influential, while attention is concentrated on 
the protagonist electrode. 

When this de-coupling operation succeeds one electrode can 
be usefully studied using the standard methods of the chemical 
reactor theory and the local electrode kinetics can be more eas- 
ily determined. On the other hand, the results of these kinds of 
approaches should often only be considered as orienting, while 
an effective and quantitative kinetic characterisation requires the 
contemporaneous consideration of both the coupled electrodes 
and their complexity in terms of velocity and composition fields, 
that is the use of complex physical—mathematical models at the 
laboratory scale. The use of detailed models is always neces- 
sary, in particular, when the thermal regulation is not completely 
effective, so that the temperature field cannot be considered 
uniform. 

Simpler cases are encountered when the experimental inves- 
tigations use perturbation methods, that is transient conditions 
following variations, often very little variations, in some opera- 
tive condition such as temperature, composition, velocity and so 
on. Some authors [12,13] have successfully used flow rate and 
composition perturbations on the electrodes of a molten carbon- 
ate fuel cell to characterise the diffusive limits of the anode and 
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cathode. This technique is now under discussion as a system- 
atic diagnostic instrument for Ansaldo Fuel Cells’ experimental 
cells. In the second part of this paper some preliminary data 
of this kind will be reported and analysed and the difficulties 
encountered in the interpretation will be underlined: in fact, 
when the variations are large enough to produce clearly mea- 
surable effects, they are often out of the linear range. Moreover, 
similar methods could be useful in the development of poly- 
meric membrane fuel cells, for instance the ENEA PEMFCs, as 
discussed in the following chapters. 

Passing from a theoretical discussion of some simple ref- 
erence cases (Part 1) to some examples of their application 
to the interpretation of experimental data (Part 2), this work 
intends to recommend a rationalisation and generalisation of 
the approach to the analysis of the diffusive resistances at the 
laboratory scale. In this way, an improvement in the diagnostic 
experimental methods can be attained for different applicative 
sectors and, in particular, more precise and reliable results in the 
characterisation of different kinds of fuel cell electrodes. 

Finally, it is to be underlined that the point of view considered 
in this work is typical of an engineering approach: as chemical 
engineers start from a known intrinsic chemical kinetics and 
work on it considering, in particular, mass, heat and momentum 
transport phenomena on both the aggregate level (e.g. solid par- 
ticles, drops, bubbles) and the reactor level, similarly a known 
intrinsic electrochemical kinetics is here assumed (e.g. the Tafel 
equation) and then the effects of transport phenomena on the 
electrochemical kinetics are discussed. As said, on the commer- 
cial scale, mass and heat transport is fully relevant and have to be 
taken into account by means of detailed physical—mathematical 
models; on the laboratory scale the attention must be focused 
mainly on mass transport phenomena, as isothermal conditions 
are granted. At both scales momentum transport phenomena 
are generally less important in fuel cell electrodes, as laminar 
regimes at high Peclet numbers are involved. Intrinsic electro- 
chemical kinetics is not discussed at all. This term, summarizing 
a number of complex molecular phenomena on a very small scale 
(transport through the ionomeric layer, activation steps, molecu- 
lar electrochemical mechanisms, and so on), is here schematized 
by means of a relatively simple relationship. On the contrary, 
emphasis is given to the question of the reference conditions for 
the electrochemical kinetics: this is a relevant point in practical 
applications because the choice of an external reference allows 
to write a more explicit and univocal kinetic expression being 
dependent on the reaction advancement. 

The usefulness in this research field of an engineering 
approach well established in chemical engineering procedures 
can be further underlined. Researchers starting from more spe- 
cific electrochemical studies could found in this similitude some 
technical suggestions, while chemical engineers could start from 
here to fully recognize electrochemical processes as integrant 
part of their field of study. The following sections would give 
some attention to this ambitious aim; also for this reason the 
discussion is maintained, especially for the Part 1, at a general 
level, without a specific reference to a particular kind of fuel 
cell or electrode and following an analytical approach to focus 
attention on these basic concepts. 


2. Local diffusive resistances 


Fuel cell electrodes are open systems in which a neutral 
reagent A, e.g. the hydrogen at the anode, must migrate from 
the gaseous phase to the reaction place. The use of Tafel’s equa- 
tion for the electrode current is standard practice, but it is well 
known that attention has to be given to the concentration effects 
of neutral reagents and to the choice of the reference conditions. 
For a “simple” electrode, defined in terms of uniform tempera- 
ture, voltage and concentrations, the choice of the reference is 
obvious; on the other hand, for a complex electrode, that is a 
“non-simple” one, at least in terms of a non-uniform composi- 
tion field, the exchange current and the equilibrium voltage are 
locally related to the composition C;, as well as to any other dif- 
ferent reference composition Cj,, for instance the composition 
Cio of the inlet gas, with particular attention to the concentration 
Cao Of the key reagent A. It is necessary, moreover, to assume 
some kind of dependence, for instance a simple proportionality, 
between the exchange current and the concentration Cas close 
to the electrode. 

The choice of a reference condition is a delicate question 
that is very often underestimated and neglected. Further details 
are reported in Appendix A. Here it should be noted that the 
choice of a unique reference composition for all the points of an 
electrode, under uniform temperature and voltage conditions, 
makes it possible to express the local current density, that is 
the local reaction rate for unit surface, merely as a function 
of the composition (or the concentration of a key reagent, or 
the utilisation factor in relation to the key). The choice of a 
unique reference condition can be extended, when needed, to 
the electrode where the temperature and/or the voltage are not 
uniform. 

It is worth noting, in particular, that the uniform voltage 
hypothesis implies a good de-coupling operation with respect 
to the other electrode: the corresponding total voltage losses are 
considered uniform or, at least, their distribution is assumed to 
be close to the mean value. On the contrary, in many instances, 
only the total voltage difference between the two electrodes can 
be assumed to be uniform, while the voltages of the electrodes 
are both non-uniform. Therefore, the de-coupling hypothesis has 
to be a priori deeply discussed and carefully a posteriori verified 
before judging the significance and the reliability of a single 
electrode analysis. 

The choice of a unique and “external” reference, if coherent, 
can be of particular aid: for instance, when the inlet gas compo- 
sition is chosen as the reference for both the reaction kinetics and 
the reaction rate and a linear and irreversible kinetic function, 
to simplify the discussion, is assumed, then the local electrode 
kinetics can be written as follows (see Appendix A) 


S 


j i P p( No) n (V — ) == (1) 
= exp (no), =(V V 
i ol Cao ° s RT 


where no expresses the electrical voltage difference between 
the electrode working under steady state conditions and a ref- 
erence electrode under equilibrium conditions at the reference 
concentration Ca. The extension to linear reversible kinetics 
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is straightforward, the extension to non-linear kinetics is only a 
little more complex. 

In Eq. (1) allowance must be made for the mass transfer resis- 
tances between the bulk fluid at concentration C4, and the fluid 
at concentration Cas close to the electrode, by expressing the 
ratio Cas/C, in terms of the limit current, that is the maximum 
diffusive flux of the limiting reagent 


Jt=keCanF, jio =keCaonF, j=k(Ca—Cas)nF (2) 


In such a way the local electrode kinetics can be written 


C : 
J = Jxo (£ ) (1 = +) exp(fno) 
Ao JL 
c : 
Siy ( a L) exp(Bio) 3) 
Ao JLo 


The second equation in (3), in particular, demonstrates that 
the local current density only depends on the local composi- 
tion. In the equations in (3) it is explicitly considered that the 
electrode can be locally in contact with a bulk fluid at a con- 
centration Ca < Cao; in turn, the concentration Ca, < Ca close 
to the electrode can be even lower due to diffusive resistances 
associated with the current. The voltage losses become negli- 
gible at very low currents and the reference concentration, but 
when the low current condition is stressed the equations in (3) 
tend to fail because they neglect the reversibility of the reaction. 
On the other hand, they are sufficiently general and reliable for 
many practical purposes and they are useful in a discussion of 
the role of transport and flow phenomena in the performance of 
the electrode. As said above, the extension to more complicated 
kinetics is only a matter of form. 

For complex electrodes working with a continuous reagent 
flow, Ca and Cas are local values depending on their position. 
In such instances, the use of the local reference Ca 

Ek ei oe a 4 
J = jx ( +) exp(B7), n = ( ea) R7 (4) 
leads to a kinetic expression that is only apparently simpler than 
the second equation in (3). In fact, in Eq. (4) the parameters j«, jL 
and 7 are all local variables, depending on their position, while 
in the second equation in (3) the local current density depends 
only on the local concentration Ca. Of course, in both cases the 
current density J referring to the whole electrode corresponds to 
an averaging operation on the electrode surface. 

However, when the operating current is much lower than 
the limit current (for instance 3j*+o <j<j_/3), the concentration 
polarisation, according to an acceptable approximation of the 
first equation in (3), can be assumed to be proportional to j/ji, 


1 j Ca J 
w n(+) in( )+- (5) 
j (3) | Jxo Cro JL 
J cost. j 
Noo ™ 7 = (6) 
© Bi Cake 
This circumstance allows a simple and direct interpretation of 


voltage—current measurements, in order to evaluate the effects 
of mass transfer resistances. In particular, when perturbation 


methods are used under constant current constraints [12], any 
ko variation can be directly detected in terms of the electrode 
voltage. However, Eqs. (5) and (6) only have local validity and 
are valid for the whole electrode only in the case of a simple 
one, while the extension to non-uniform composition electrodes 
requires an averaging procedure. 

Moreover, (5) demonstrates that the local voltage loss is 
partly imputable to diffusive effects (Cas < Ca) and partly to 
a non-uniform composition field (Ca < Cao): the measurement 
of No and its variations cannot be directly connected to the dif- 
fusive resistances unless there is information available to allow 
appropriate corrections to be made to the results in terms of 
non-uniform concentration fields. 


3. The analogy with heterogeneous chemical kinetics 


If, for simplicity, the discussion is limited to irreversible linear 
kinetics and the following definitions 
Jo j 

ko = ex , r= J 

o (=) p(Bno) AF (7) 

are used, so that an “intrinsic” kinetic constant k, of the elec- 

trochemical reaction is defined as a function of the electrode 

over-potential, the classical formulation of a first order hetero- 


geneous reaction on the boundary is obtained instead of Eqs. (1) 
and (3): 


r = koCas = ke(Ca — Cas) 


= KC K= ae k (8) 
= KCa, Ae ae 


In local terms, or in terms of a simple electrode, the apparent 
kinetic constant and the limit kinetic constant coincide with the 
overall constant and the transport coefficient, respectively 


(a-2% kas k (9) 
and the relationship between the apparent (ka), the limit (kg) 
and the intrinsic (ko) constants, equivalent to Eq. (4), is the well 
known 


ko/ kL ) (10) 


Ka ( 
ky \(+ko/kt) 


or, in the form more common in electrochemistry, 


ae oi 
= (£2) expire) (1 +) d1) 


For the simple electrode, according to Eq. (11), the depen- 
dence of the intrinsic constant ky on the over-potential can 
be obtained from low current measurements (3J*9 < J < JL/3), 
where Cas = Ca and ko < ky (negligible diffusive resistances); 
in turn, the limit constant can be obtained from high current 
measurements, where ko > ky (negligible reaction resistances). 
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Eqs. (11) and (10) are different obviously only for the form, 
but they correspond to two different points of view, which must 
be reconciled in order to appreciate some other fundamental 
analogies. In heterogeneous chemical kinetics the constant ko 
is considered as a real constant at constant temperature, it is 
often considered as not being directly measurable but deducible 
from global kinetic measurements (ka) and transport coefficient 
evaluations (kc). A variable and directly measurable intrinsic 
kinetic constant should be automatically felt as strange, even 
if, in principle, it could correspond to controlled temperature 
changes at the reaction site. On the contrary, in the electrode 
analysis, the control of the intrinsic kinetic constant can easily be 
obtained through voltage measurements: while ka is still directly 
connected to global kinetics, that is the electrode current, ko is a 
variable directly connected to the electrode over-potential by a 
kinetic equation like the first one in (7) and the measurement of 
both allows an evaluation and an analysis of the mass transfer 
phenomena occurring at the electrode in terms of the constant 
k, that is the limit current. 

After clarifying the different points of view, the analogies 
existing between some electrochemical kinetic problems and 
some classical heterogeneous chemical kinetic formulations can 
be reviewed, so providing wide and well-established literature 
references. In particular, the recourse to such analogies in the 
averaging operations of complex electrodes often provides the 
opportunity for utilising in this case solutions already known in 
other contexts. A well known classical example is the analogy 
between porous electrodes and porous chemical catalysts [14]; 
here another analogy will be specifically used, the one between 
the kinetic behaviour of a planar electrode of finite dimensions 
in contact with a moving fluid and the well-known problem of 
diffusion with a chemical reaction at the boundary. 

The relevance in terms of applications is confirmed taking 
in account that, even if fuel cells usually work at not too high 
current density, also under apparently safe conditions a large 
part of the cell can locally work in conditions where diffusion 
phenomena significantly affect performance [11]. 

Other analogous approaches that refer to the chemical reactor 
theory can be used when the optimisation of the electrode shape 
and the flow pattern in an electrochemical reactor are considered. 
As examples, the studies of the comparison between interdigi- 
tated and serpentine configurations in PEMFC electrodes can be 
cited [15,16], a topic which seems to be worth further discus- 
sion elsewhere. Here the comparison between a well-mixed and 
a longitudinal flow electrode will be discussed in some detail 
(see Part 2, Appendix D) and the results of the longitudinal-flow 
electrode study will be used in the analysis of MCFC anodes 
(see Part 2). 


4. From local to global kinetics 


Eq. (11) is valid in local terms, but it is usually extended to 
much wider contexts. It is worth exemplifying the conditions 
under which it can be correctly extrapolated. 

For the whole electrode, the definition of the apparent and 
limit constants has to be made in terms of the mean values of 
the reaction rate (or current density), obviously with reference 


to the inlet concentration: 


k= m J 7 KCam 
2 Cao nFCao Cao’ 
J kC 
kL TmL L _ Ke AmL (12) 


~ Cao i? nFCho Cao 


As a consequence, the local Eq. (11) can be quite inadequate 
and require a correction. By defining 


1 Cam UfpCads) 


= = : >1 (13) 
A Cao (CaoS) 
the preceding Eq. (12) can be summarised, for linear kinetics as 
Ss ee. (14) 
a— A $ G= AL 


It is useful to compare a complex electrode, working at a non- 
uniform composition, with a corresponding simple electrode, 
working at uniform concentration C4o; for the simple electrode 
the local Eq. (11) can be written as 


1 1 1 1 


1 1 
S— + = - = | — J exp(— no) + — (15) 
K ko ke J (=) i j JL 


For a complex electrode, instead of the local constant K, con- 
taining the reaction constant ky and the transport constant ke, 
reference can be made to an averaging operation on the electrode 
surface 


1 1 1 1 1 


Z 16 

Aka Ko T ke ko T ALkL ao 

: DE ea a7) 
= ex ——— 
AJ Jxo P ALJL 


so that the apparent constant ką and the limit constant Ay can 
be connected to local constants through the introduction of the 
coefficients A and AL, which take account of all the effects of a 
non-uniform concentration field. 

The validity range of Eq. (16) is wider than that of Eq. (15), 
but it has its limits. In the first place, as said above, ky must be uni- 
form, that is the electrode must be isothermal and efficaciously 
de-coupled. In the second place kę must also be uniform: partic- 
ular shapes or flow fields can be associated with local differences 
in the transport coefficient; this, moreover, can be affected by 
local non-linear phenomena such as those connected with the 
bulk flow towards the electrode (see Part 2, Appendix E). 

However, in many instances, Eq. (16) can be taken as a guide 
for a better understanding of some typical operations in the 
analysis of chemical and electrochemical reactions. In heteroge- 
neous chemical reactions reference can be made to the following 
typical identification problems: 


e On the local scale, that is when it is possible to work on sim- 
ple and uniform systems, the characterisation of the intrinsic 
kinetics (kinetic constants and reaction orders) is relatively 
easy. 

e Ona larger scale, that is in terms of aggregates or reactors, the 
characterisation of the intrinsic kinetics passes through a more 
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articulated physical-mathematical description of the system: 
the global or mean kinetics (ką) are measured, the transport 
coefficients (kc) are estimated, the composition fields (A) are 
simulated and averaged and the intrinsic kinetics (ko) are, 
finally, determined. Only in rare cases the intrinsic kinetics 
are already known from preceding measurements on simpler 
systems, so that the transport phenomena can be the specific 
objective of the identification problem. 


When dealing with electrochemical kinetics, the terms of the 
discussion remain substantially the same only when the local 
scale is under consideration (laboratory systems which can be 
considered simple electrodes), but they become rather different 
at greater scales (complex electrodes), where the identification 
problem involves the measurement of the mean kinetics (J, ka) 
and intrinsic kinetics (V, ko) as well as the interpretation of the 
limit kinetics (JL, kL) in terms of transport coefficients (kc) by 
taking into account the concentration distribution (A, AL). The 
identification problem is often complicated by other unknowns 
(symmetry factor 6, reaction orders) appearing in the relation- 
ship between V and ko. 

In both cases the preliminary measurements are restricted, 
for what is possible, by well-controlled laboratory conditions, 
where the identification problem can be simplified by reference 
to uniform concentration fields (A = Ay = 1) and non-controlling 
transport phenomena (ke >> ko). In both cases the experimenta- 
tion, at the laboratory or pilot scale, is then extended to systems 
of larger sizes and greater complexity and, finally, the identifi- 
cation problems give way to prevision problems, in particular 
design problems relating to developmental and commercial pur- 
poses, on systems where the approximations under discussion 
are nearly always unsuitable. 

In other words, and to stress the point, Eq. (16), or, in more 
electrochemical language, Eq. (17) shows that the identification 
of the global transport coefficient of an electrode, ke, requires 
a number of operations that are not always straightforward or 
certain and, in any case, a non-superficial level of theoretic 
analysis. 


e Inthe first place, the identification is possible by means of the 
consideration of only one electrode on condition that the other 
electrode is working under an approximately uniform voltage 
field. When this is not the case, the consideration of only one 
electrode cannot be considered significant and the transport 
phenomena of the two electrodes cannot be de-coupled. 

e In the second place, the experimental work can be done 
according to essentially known methods, by achieving cou- 
ples of values for the current (J, that is ka) and the voltage (no, 
that is, substantially, ko). 

e Finally, an estimation of the mean electrode concentration 
(Cam, that is the parameter A) has to be made for the various 
experimental conditions used and, in particular, for the limit 
current conditions (JL, AL). The evaluation of these means 
can be done experimentally, from local concentration mea- 
surements taken at various points on the electrode, but an 
alternative to this experimental burden can be found in the 
calculation, starting with few known concentration values: 


usually, the inlet value Cao and the outlet value (1—w)Cao 
are at least known, but employing a simple arithmetical mean, 
such as Cam = Cao(1—w/2), may be inadequate. 


In addition, particular attention has to be given to the depen- 
dence of the transport coefficient on the flow conditions of the 
reagent fluid in relationship to the electrode. The effects of the 
flow conditions are usually appreciable on all the measured or 
estimated parameters, especially the transport coefficients and 
the averaging coefficient A. In the simplest cases the flow con- 
ditions can be taken into account in terms of a further kinetic 
parameter alone, the “flow dynamic” constant 

q 
kf = 5 (18) 
that is the fluid flow rate per unit electrode surface, to be put 
beside the previously considered reaction and transport con- 
stants (see Appendix C1). 


5. Analysis and correlation methods 


The approach discussed here highlights the advantages of a 
unique and external reference. This approach will be underlined 
in Appendix A. 

Some considerations on dealing with the organisation and 
the correlation of the experimental results according to what has 
been reported in the previous section are presented in Appendix 
B. In Appendix C some examples of a theoretical nature are 
discussed, with the aim of illustrating the real relevance of the 
averaging operations. 

In Part 2 of this work different analysis methods will be 
discussed in relationship with different experimentation tech- 
niques. In addition, concrete examples of their application 
to Ansaldo Fuel Cells data will be presented in detail. In 
Appendixes D and E of Part 2, moreover, the effects of the flow 
field on the electrode performance will be discussed on the basis 
of the chemical reaction theory and its classical limit models. 
Many more examples could be added, but the ones cited can 
be considered sufficient to provide a comprehensive overview 
and to draw some conclusions and rules for a more exhaustive 
and reliable analysis of the experimental current-voltage data in 
terms of diffusive resistances. 


6. Conclusions 


The analysis of the diffusive resistance is made difficult by 
the non-linearity and complexity of the phenomena even at the 
laboratory scale. 

The problem can be tackled at the experimental level by using 
a number of different techniques: differential or finite perturba- 
tion methods, constant flow or constant current methods, and so 
on. However, the effects of each variation cannot be interpreted 
as a simple phenomenon, but the result of at least two factors: 


e the local diffusive resistances set and 
e the concentration distribution of the electrode. 
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The first can only be determined after estimating the second 
by means of an adequate physical—-mathematical description. To 
achieve this it is particularly useful to make recourse to a close 
analogy between electrode kinetics and heterogeneous chemical 
kinetics with the reaction at the boundary and to use this analogy 
within the chemical reactor theory. 

As Part 2 will show in detail, this kind of analysis can be 
successfully applied to fuel cell experimental data. 


Appendix A. Kinetic expressions and reference 
conditions 


Given the local electrochemical kinetics of the kind 


j = KITCK exp(Bos) — K'TCiy expl(B — Dvs], 
VnF 
v= — 
RT 

which contains the concentrations and the voltage acting on an 


electrode point (Cis, vs) and must satisfy the thermodynamic 
consistence conditions 


k 


k! 


(A.1) 


iC?! exp(—veq) = MiC exp(—ves) 
= TI;C}i exp(— Ver) (A.2) 


a reference can be chosen at the equilibrium condition corre- 
sponding to an arbitrary composition Cir 


Jer = KIT;C% exp(Bver) = KTC® expl(B— Dvel (A3) 
so that 
oo, Cis \“! are aed 
J Jarl i ( c) exp(fnr) i Ca 
exp[(8 — 1)n:l}, r = Vs — Ver (A.4) 


In such a way, the local kinetics is written in terms of the 
voltage loss in relationship to the equilibrium voltage of the 
chosen external reference condition. This advantage has a price: 
the local voltage loss in respect to the bulk fluid at composition 
Ci, with which the irreversible phenomena of energy dissipation 
are associated, is different from the one appearing in the kinetic 
expression (A.4) 


N = Vs — Veg = Vs — Ver + Ver — Veq 
Cir 
= Nr + Xvi In ($) 


As the electrode concentration Cis and the local concentration 
of the bulk fluid C; are connected in transport equations, for 
instance under the form 


(A.5) 


1 nF vij 
j = —keci(Ci — Cis) 7 Cis = Ci + 


A.6 
i n Fkei ( ) 


in Eq. (A.4) the concentration difference due to mass transfer 
can be taken into account in terms of the transport coefficient kg; 


or in terms of limit currents, with reference to local conditions 
(Li) or to reference conditions (jr) 


; nF . nF 
jti = —kciCi—, Jiri = —kei Cir — (A.7) 
Vi Vi 
yielding 
Cis _ Ci (1+ vij )=£4 vij _ Gi d 
Cir Cir n FkciCi Cir  nFkca&Cir Cir Jiri 
(A.8) 
and 
ae Ci vj 1" 
S a de 
J Jer Et] 
ailhenn oy 
Paea Cy T e 


Besides the reaction order (a;) and the other coefficients of 
a stoichiometric nature (v;, §), this final form of the kinetic 
expression (A.9) contains 


e the reference composition Cir; 

e the exchange current j»; at the reference composition; 

e the transport coefficient between the fluid and electrode kei 
or, alternatively, the limit current jur; at the reference compo- 
sition; 

e the voltage loss 7; in relationship to the reference equilibrium; 

e the local composition Cj. 


In particular, in many instances, the variable n, can be con- 
sidered uniform on the electrode while, on the contrary, the local 
variable 7 surely is not. Moreover, the non-uniform concen- 
tration variable C;, in many cases of practical interest, can be 
expressed through the definition of a unique degree of advance- 
ment 


Ci = Cir + vix) (A10) 

or the utilization factor of a key reagent A 

Ca=Ca(l—u), C= Cr— (=) Caru (A.11) 
VA 


Therefore, Eq. (A.9) is equivalent to a kinetic function, at con- 
stant temperature and voltage, depending on an advancement 
variable alone (x or u). 

An electrode can be defined “simple” if all the local variables 
(T, P, nr, x) are independent of their position, that is each assumes 
a unique uniform value on the electrode. For many complex 
electrodes, especially at the laboratory scale, the temperature, 
the pressure and the voltage can be still assumed to be uniform, 
while the composition (degree of advancement) cannot. 

It is worthwhile commenting further on the electrode voltage. 
At uniform temperature and pressure the two coupled elec- 
trodes have, at least at the laboratory scale and with a reasonable 
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approximation, equal current and a uniform voltage difference: 


fi x)dS = if I (np, x) dS", 
s s! 


(JS = J'S), N: + ni. = cost (A.12) 


When one of the electrodes is considered to be simple (for 
instance the one marked with the apex), the uniformity of 1}. 
implies the uniformity of 7; on the other electrode, although it 
is not considered simple. In such cases, electrode kinetics of the 


type 
j= FH) 


allow the extension of many of the identification, design and 
optimisation procedures of the classical theory of isothermal 
chemical reactors to electrochemical reactors. 

According to Eq. (A.9), a particularly simple and interesting 
case to analyse as a reference example is the one regarding a 
unique limiting reagent A, with reaction order and stoichiomet- 
ric coefficient equal to unity (wa = —va = 1). Then the kinetic 
expression becomes 


(A.13) 


| Jer exp(Bnr) J 

= C Car — A.14 
J | Ca | A nF ke ‘Ar €Xp(—nr) ( ) 
and, if 
Cae = Car €Xp(—Mr) (A.15) 


is the concentration at which the reagent fluid is in equilibrium 
under a given over-potential nr, the electrochemical kinetics can 
then be reduced to simple linear kinetics 


.__ | Jar exp(Brr) J 
= C C A.16 
J | Ca | a — Cae — TFK (A.16) 
By putting 
J r Jur EXp(Bnr) 
Pe K = ’ o> 
nF (Ca — Cre) nFCar 
(A.17) 
Eq. (A.15) becomes simply 
re a 
r=K(Ca— Ca),  K=(—+— (A.18) 
ko ke 


A similar result can also be achieved when more complex 
kinetics are considered. For instance, if two or more reagents 
and non-linear dependences are involved, an expression just like 
the preceding is derived in the range of validity of a linearisation 
of Eq. (A.9) for low utilisation factors (u< 1) and low currents 
G <j Lri ). 


Appendix B. How to determine the complexity 


When the intrinsic kinetic constant kọ increases, the appar- 
ent kinetic constant ką is limited by the limit constant k, 
(0 < ka < kL). Therefore, both variables 


k 
X=, 


i (0< X < 1) 


(B.1) 


yo e- (0<Y<1) 


5 (B.2) 


can vary in the range between zero (kj =k,=0) and unity 
(ko >> kL; ka =k). 

In local terms, or for a simple electrode with linear kinet- 
ics, the limit constant coincides with the transport coefficient 
(kL = kc) and the apparent constant depends only on the intrinsic 
constant and the transport coefficient (1/ka = 1/K = 1/ko + I/kc), 
so that between X and Y the simple relationship 
Y=xX (B.3) 
exists, that is equivalent to Eq. (15) and in an XY plot this 
corresponds to the square diagonal. 

On the contrary, for a complex electrode, with linear kinetics, 
the composition distribution on the electrode, with Ca <Cao, 
implies that the apparent constant as well as the limit con- 
stant is less than the respective reference values (kL < kc; ka < K). 
Therefore, Eq. (B.3) must be correct, for instance in the form 


Yo _ (ke — kt) 
x7 = f(X), fe= 


ia ee: WA 
k +( ) ~ (B.4) 
that is equivalent to Eq. (16). 

When the corrective function fe is null for all the values of 
ko, Eqs. (15) or (B.3) are recovered, while Eq. (B.4) with fe >0 
corresponds to a curve which is under the diagonal as much the 
non-uniformity of the diffusive phenomena taking place at the 
electrode is important. 

Of greater interest are the asymptotic solutions to which Eqs. 
(16) and (B.4) degenerate at the extremes of the range. At one 
extreme, in the field of low currents, a proportionality between 
X and Y can be assumed 


Y 
ko & kL, Y=coX, co=lim š] (B.5) 
X>0 


X 

and the value of the constant co can be obtained as an asymptote 
of the experimental XY plot or in terms of the derivatives of ką in 
relationship to ko. If X= Y=0, ko/ka = 1 and dka/dko also = 1. In 
the neighbourhood of this point the difference between ka and 
ko is only appreciable in terms of the second order derivative, 
thatis ka = ko + (d?ka/dk2) k2/2. Therefore the proportionality 
constant is 


(B.6) 


At the other extreme, in the high current field, the propor- 
tionality can be assumed to be between 1—X and 1—Y 


ko > ke, 1— Y = c(l — X), 


: fed 
cL = lim 
X>1 


ess (B.7) 
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and the proportionality constant s 
ki 
CL = Fan, alt 
(ko dka/dko) x,t, 


—1 
= [Ar + (1/kXkdA/dko u] © far =(L- D 


(B.8) 


as k, becomes indistinguishable from k, for high, even if finite, 
values of ko, the limit of kždka /dko for ka = k is often zero and 
c diverges, that is the function Y(X) tends to the vertical X= 1. 

Another parameter which can generally express the differ- 
ences between Eqs. (B.3) and (B.4) is the size of the belly of 
the curve, as measured by the segment on the other diagonal 
(Y=1 — X; ko = kL) of the square, between the diagonal crossing 
and the intersection with the curve. The length dy of this segment 
is 


_ (42/2) fo 

(2— fo) 
where f.# indicates the correction in Eq. (B.4) for ko =k. 
Obviously, when Eq. (B.3), that is Eq. (11), is correct, the con- 


stants assume the reference values corresponding to the principal 
diagonal Eq. (B.3): 


dy (B.9) 


Co=l, cL = 1, dz =0 (B.10) 


whereas when there is some shifting from the diagonal, the 
constants will, in most cases, be Eqs. (16) and (B.4): 


Co < 1, cL > 1, dy > 0 (B.11) 


The low current constant co is particularly important for the 
proper characterisation of an electrode, as it directly modifies 
the linear current—voltage relationships of type (5). 


Appendix C. Some examples 


The real relevance of the asymptotic constants cy and cy can 
be illustrated by reporting some significant examples. 


C.1. Stirred electrode 


A first, very simple, example regards the behaviour of one 
well stirred electrode to which the reagent fluid is continuously 
fed, so that a unique value of concentration Ca is sufficient 
to describe its steady condition (open system with concentrate 
parameters, or continuous stirred tank reactor, CSTR). So, the 
local condition is described by Cag, while the inlet concentration 
is Cao, which is connected to the former by the balance equation 

K 
Cao-Ca=—Ca, ktsg (C.1) 
kg 
and the flow dynamic constant ks, (volumetric flow rate supplied 
to the electrode per unit active surface). 
In terms of the electrode concentration, according to Eq. (C.1) 


Cao 


CA = G4 Kk 


(C.2) 


the apparent constant referring to the inlet concentration 
ka = KCa/CaAo is then 


ka = = - (+i) (C.3) 
° (K/K) \ko ke ke 
and corresponds to an averaging coefficient 
K 
A=1+— (C.4) 
kg 


Eq. (C.3) seems different from Eqs. (10) and (11), at least 
because of the flow constant kp. In fact, if the discontinuous 
concentration drop from inlet Cao to internal Ca is taken into 
account, the electrode should be classified as “complex” and 
behavioural differences are to be expected when it is compared 
with the simple electrode at uniform C4o: the X—Y trends of a 
stirred electrode with the position kL = ke typical for the simple 
electrode should be represented by a family of straight lines with 
increasing slope with increasing kf. Nevertheless, from another 
point of view, the stirred electrode can be considered “simple” 
because all its surface points work at uniform concentration CA. 

In order to reconcile the two antithetic positions it is sufficient 
to observe that, in the absence of reaction resistances, the limit 
constant is 


tal eee C3) 
ca? ae ' 
which is equivalent to 
ke 
Ata (C.6) 
kf 


and also depends on kf; in such a way it is possible to re-obtain 
Eqs. (10) and (11) through the elimination of ke in Eqs. (C.3) 
and (C.5). 

Eqs. (C.3) and (C.5) come from definition (12). In other 
words, the validity of Eq. (11) is extended thanks to defi- 
nition (12) of the apparent and limit constants, a definition 
that has been stated in direct agreement with the experimental 
evidence. 

However, it is worth noting that the simple identity between 
the transport coefficient and the limit constant in the second 
equation in (9) is not adequate any more, while the valid position 
has become the more complicated Eq. (C.5), in which the “flow” 
constant kf is also considered in order to take into account the 
limits imposed on the electrode current by the amount of reagent 
at its disposal. Obviously, only if the link between ky, and ke is 
known, for instance the one in (C.5), Eq. (10) or (11) makes it 
possible to analyse the experimental results on the limit current 
(kL) in order to achieve a reliable characterisation of the transport 
phenomena (ke). 


C.2. Longitudinal flow electrode 


A second example, which is apparently just as simple as the 
preceding, but able to highlight the possible difficulties of the 
problem, regards a planar electrode in contact with a longitudi- 
nally flowing and transversally well mixed fluid: it is the system 
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characterised by a piston flow, that is the one defined as an ideal 
tubular reactor (or plug flow reactor, PFR) in the chemical reac- 
tor theory. Under stationary and constant flow rate conditions, 
the differential mass balance along the axial coordinate z, written 
in terms of the utilisation factor 


kedu’ = K(1 — u')dé, ¢=0 w=0, ¢=1 w=u 
1 1 

a= =T tat 

allows the integral 

K = —k¢ ln(1 — u) (C.8) 

and the limit condition 

ke = —ke ln(1 — uL) (C.9) 

Moreover, it is 

keu = ka, ke = uLkL (C.10) 


and, by combining Eqs. (C.8)-(C.10) with the last equation in 
(C.7) the variable X and Y are obtained in terms of the variable 
u used as a parameter 


ka u 
X= m = — 
sa ús i (C.11) 
fol] =i X 
p TAA] pa M-m 


This time the similarity between the second equation in 
(C.11) and Eq. (11) is only very partial. Obviously both equa- 
tions, thanks to the way they have been constructed, allow the 
limits where the reaction (ka =0 for ko < kL) or the transport 
(ka =k for ko >> k) is controlling, but the characteristic con- 
stants of the curve (B.4) 


_ (kv) (ke + 2ko) 
2k. 


_ (kL) [ke 
a= (5) laa 


eo 
kL ko=ki 


(C.12) 


depend on the ratio kc/k¢ (see Figs. C.1 and C.2) and assume the 
values (B.10) only in the trivial case in which the flow constant 
is quite non-limiting, or non-influencing (kf >> kL). Under such 
conditions the amount of the fed reagent is much more than 
that consumed (very low utilisation factors) and the reagent 
concentration is almost uniform along the length of the elec- 
trode, so that the whole electrode behaves as a simple electrode 
(Ca = Cao). 

Under the contrary conditions, when the amount of the 
inlet reagent becomes controlling (kf < kL), the characteristic 
constants (C.12) approach the limit values co = 1/2, 1/cĻ =0, 
respectively, and Eq. (C.11) tends to degenerate into ka/kL =1, 
which is quite different from Eq. (11). 

A first comment, of a reassuring nature, could be that a 
deliberate choice of flow-controlling conditions (kf < ke) is very 
unlikely, so that the values assumed by the asymptotic constants 


in (C.12) are almost always close to the reference ones in (B.10): 
in such instances Eq. (11) is substantially acceptable for all the 
values of ko, that is the voltage, and the electrode behaves as a 
simple one. 


) -1 (C.7) 


Just this kind of consideration, on the other hand, diverts 
attention from the circumstances in which the Eqs. (C.12) are 
more pertinent, so impeding a full appreciation of their role in 
possible malfunctions, often not foreseen nor suspected, in a par- 
ticular region of a cell or stack, at the commercial or laboratory 
scale. In fact, even if usually the operative conditions involve a 
current decidedly lower than the limit one in terms of average 
values, some section of a cell can really work under conditions 
very close to the limit current and the complete exhaustion of 
the reagents. 


(0 5 10 15 20 25 30 


Fig. C.1. Longitudinal flow electrode. The asymptotic constants co and 1/c; as 
a function of the rate ke/kt. 
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Fig. C.2. Longitudinal flow electrode. XY plot for different values of the ratio 
k./k¢. The diagonal represents the simple electrode. 
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Moreover, and it is perhaps the most important point, the first 
equation in (C.12) stresses the necessity of using a corrective 
Co constant different from unity in the interpretation of many 
experimental results, even if they are obtained under low current 
conditions. 


C.3. Diffusive electrode under uniform velocity 


Another example, which is particularly relevant for polymeric 
membrane fuel cells in their special interdigitated or serpentine 
flow configurations [15], is the that of a planar electrode in con- 
tact with a fluid flowing over its surface with a uniform velocity 
field. This uniformity is typically due to a porous medium (Darcy 
flow) and the reagent can reach the electrode exclusively through 
a diffusive path orthogonal to the electrode direction. If v is the 
velocity, s the thickness of the porous medium and L its length 
and reference is made to the pellicular regime where the resi- 
dence times L/v are much lower than the diffusion times s2/Da, 
the mean electrode current can be evaluated by integrating the 
partial derivative system 


aCA a°Ca 

— = Da — C.13 
az j 3y? (C13) 
z= 0; Ca = Cao (C.14) 

ac 
y=0, Se = koCa (C.15) 

dy 

DAL 1/2 

yo ow, — Ky<s|, Ca = Cho (C.16) 


In particular, the last boundary condition corresponds to the ref- 
erence to a semi-infinite medium, which is fully coherent with 
the pellicular regime assumption. 

The operations of integration and calculation of the mean flux 
towards the electrode, after having defined the flow dynamic and 
transport constant as 


k= k= =A, (= « 1) (C.17) 
yield [16,17], for the apparent constant, 
L i) 1/2 
T 
Brk fı exp ( ks ) et Fecal (C.18) 
ko keke (kek,)!/? 


When ko/(kek¢)!? >> 1 (transport controlling) the preceding Eq. 
(C.18) can be written 


Akgk. \ 1/2 
kL = (=) 
TT 


so that, by eliminating k, between Eqs. (C.18) and (C.19) 


tk? 4k2 4k2\ 
ka = kL 1 — exp erfc 
4ko rk? rk? 


(C.20) 


(C.19) 


The characteristic asymptotic constant can so be obtained 
from Eq. (C.20) so that 


where the asymptotic solution with controlling kinetics 
(Ko/ke & 1) requires a series expansion of both the exponen- 
tial and the error functions, while for the transport controlling 
solution (koọ/ke > 1) the well known approximation expa?) 
erfe(x) © 1(/xx) for x > 1 is sufficient. 

It is worth noting here that a non-controlling feed rate con- 
dition is encountered (k,/kp << 1), where the reference values 
(B.10) should be expected to be correct, while Eq. (C.21) and 
Fig. C.3 show the necessity of a certain, non-negligible cor- 
rection. In particular, the first equation in (C.21) indicates that 
the simple electrode solution needs a correction by means of a 
coefficient co less than unity even in the low current field. 

Naturally, Eq. (C.18), as well as the other derived Eqs. 
(C.19)-(C.21), fail when the feed-rate limitations begin to act, 
that is the ratio k,/k increases until it becomes comparable with 
or even greater than unity. The complete solution for all the 
values of ke/kp = Da L/vs? can be found in specialised books 
[17]; here it is sufficient to note that for k./kg>> 1 the semi- 
infinite layer hypothesis is no longer realistic and the diffusion 
acts throughout all the layer thickness, so that other solutions 
become sufficiently accurate instead of the one presented in this 
section, for instance those discussed in [16] or that of the pre- 
ceding section for a transversally well mixed flow. For instance, 
Eq. (C.18), just as in the preceding example, also tends to finally 
be reduced to 


ke/ke > 1, 


Bek.  ka=kL=kf (C.22) 
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Fig. C.3. Uniform velocity electrode. XY plot for ke/ke & 1. 
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Fig. C.4. Porous electrode, local behaviour. XY plot for SpD pett! SexDaex = 3. 


C.4. A local example 


The examples in subsections C1 and C2 refer to effects at 
the electrode scale, for instance a “postcard” laboratory elec- 
trode. In subsection C3 attention is shifted to a level which is 
more properly classifiable as an aggregate scale: for instance, 
the considered flow pattern interests a small porous region, 
1 or 2mm in length and one or two tenths of a millimetre 
in thickness, between two adduction channels of an interdig- 
itated PEMFC. However, in all these examples, the diffusion 
and reaction paths are considered to take place strictly in 
series. 

On the other hand, from a local point of view, the porous 
structure of an electrode may act by combining diffusion and 
reaction in a series-parallel mechanism. Here the scale consid- 
ered is the one of the thickness of an electro-catalytic layer, that 
is, typically, less than 100 wm. 

In such instances, a rough and non-univocal local description 
can do reference to an external diffusive layer, corresponding to 
the external transport coefficient 


D Aex 


(C.23) 


Keex = 
Sex 


and to an internal reaction affected by diffusion, which can be 
treated in terms of the effectiveness factor g and the Thiele 


module @ 


f tgh® ka \'? 
koeff = kop = kyapSpY, Y= ee P = Sp Dee 
e 


(C.24) 
The apparent constant is then 
: : F : k k (C.25) 
ka kop kt > L = Kcex ` 


and, between the low current range, where kinetics is controlling 
(ka = ko), and the high current range, where the external transport 
is controlling (ka = kL), there may be an intermediate range where 
both internal diffusion and reaction act through the geometrical 
mean ka = (koDaetp/s)". 

Moreover, due to the series-parallel mechanism, the sum of 
the diffusive resistances decreases with increasing current, so 
that the X-Y plot lies over the square diagonal (see Fig. C.4). 

Such effects will combine with the effects acting on the other 
scales, so that an exhaustive separation and identification of the 
various effects could prove to be extremely difficult. 
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